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Abstract. Recently (see quant-ph/0503040 I an explícit example has been given of a 'PT-symmetric non- 
diagonalizable Hamiltonian. In this paper we show that such Hamiltonians appear as supersymmetric (SUSY) 
partners of Hermitian (hence diagonalizable) Hamiltonians and they can be turned back to diagonalizable forms 
^ ■ by appropriate SUSY transformations. 
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1. It is well-known that there exist non-Hermitian Hamiltonians which cannot be reduced to 
a diagonal form by the change of the basis (so called non-diagonalizable Hamiltonians, see e.g. 
[ï]). To illustrate better our ideas we will consider here only regular Sturm-Liouville problems. 
The set of eigenfunctions of a non-diagonalizable Hamiltonian is not complete in corresponding 
Hilbert space (21 El- The characteristic determinant has múltiple roots. Together with any 
eigenfunction with a simple eigenvalue coinciding with a múltiple root of the characteristic 
determinant there exists a set of associated functions [21 Ej • The linear hull of the eigenfunction 
and corresponding set of the associated functions forms for the given value of the energy the 
root subspace (see e.g. |3]). Recently an explicit example of an exactly solvable "PT-symmetric 
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non-diagonalizable Hamiltonian was given jl]. 

We have discovered that supersymmetry (SUSY) transformations may convert an Hermitian 
(hence digonalizable) Hamiltonian to a non-diagonalizable Hamiltonian, which in particular 
can possess the VT symmetry, and vice versa. The possibility which does not appear in the 
linear àlgebra. It is related with the possibility to "create" more than one "bound state" at a 
given non-degenerate value of the energy. Since the energy level is non-degenerate the other 
state cannot be an eigenfunction of the Hamiltonian but it can be an associated function. In 
the opposite process, when we are "deleting" an eigenfunction having a non-zero associated 
function, the latter is transformed to a "real eigenfunction". This looks like it "emerges from 
the background" and, therefore it may be called "background eigenfunction" . 

2. We have found that the possibility described above appears if second order SUSY 
transformations or higher are used. 

Let us consider two ordinary second order differential equations 

(h - E)ip E (x) = h = -dl + V {x) x E [a, b] (1) 
(h 1 -E)tp E (x) = h 1 = -dl + V 1 (x) xe[a,b] (2) 
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with a and b being finite numbers. We say that the Hamiltonian hi is related with h by a 
second order SUSY transformation if 

1. V l = V -2[logW(u 1 ,u 2 )]" (3) 
W(ui, u 2 ) = uiu' 2 — u[u 2 7^ Vï 6 (a, b) 

2. (h - a l>2 )u lj2 {x) = 0. (4) 
It is known (see e.g. jHj) that in this case tp E is related with ipE as follows: 

íp E = Lip E = W(u 1 ,u 2 , ^b)/W(ui, u 2 ) E^a u a 2 (5) 
Voa.a = «a.i/W'^i, w 2 ) • (6) 

We are using the symbol W to denote Wronskians and will everywhere suppose that W(ui, u 2 ) ^ 
Vi e (a, b). Equation © holds for any ip E from the two dimensional space ker(/i — E). The 
operator L intertwines the Hamiltonians h and hi, Lh = h\L. 

Let us suppose that Vq(x) is a real-valued and sufficiently smooth function for x G [a,b], 
Consider two boundary value problems, that we will denote (I) and (II) respectively, defined 
by the equations (JTJ) and (J2J) and the boundary conditions 

Ma) = Mb) = (7) 

¥»s(o) = ^(6) = 0. (8) 

It is well-known (see e.g. jïïj) that the problem (I) has only discrete, simple and real spectrum 
of eigenvalues E = E n , n = 0, 1, 2, . . .. 

We will now formulate conditions for u\ and u 2 leading to a complex- valued V\(x) given by 
with a real and simple spectrum coinciding with the spectrum of Vo except for one level 
and the Hamiltonian h\ is non-diagonalizable. 

It follows from (J3j) and (0) that 

Ve = T777- \ [( w ( u ii U 2)E + a 2 u[u 2 - aiUiu' 2 )ip E + («i - a^Mi^í?] ■ ( 9 ) 

W(ui,u 2 ) 

It is clear from here that if both u\ and i/je, E ^ ati,a 2 , satisfy the boundary conditions (JTj) 
then íçe given by (Jüj) satisfies the boundary conditions (JHJ). The only possibility for u\ to satisfy 
the zero boundary conditions is to be an eigenfunction of h , u\ = ipE k , so that it is (up to 
an inessential phase factor) real and a\ = Ek, which we shall suppose to be the case. This 
means that the Hamiltonian h\ has the same spectrum as h except maybe for the vàlues ot\ 
and a 2 but since u\ is supposed to satisfy the boundary conditions (J7J), the function ip a , 2 given 
in (jüj) is an eigenfunction of h\ and E = a 2 is the spectral point for h\. Remembering that 
we want to keep the real character of the spectrum of h\ we have to choose a 2 real also. So, 
we choose both a± and a 2 to be real and the function u\ is fixed to be real but we want to 
get a complex potential difference defined by equation (jïïj). This is possible if u 2 is a complex 
linear combination of two real linearly independent solutions of equation (JTJ). Let a 2 (^ ai) 
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also coincides with a spectral point E\ of ho, a 2 = E h and u 2 = i^Ei + icifj E , ceR where ipE l 

satisfies the boundary conditions ((Zj) and ip J is any real solution of Eq. ((TJ) at E = E\ linearly 
independent with ipE r We notice that u 2 (x) ^ Va; G [a, 6]. 

We claim that with Ui and u 2 being chosen as it is described above the potential V\ given 
in (jïïj) has the spectrum coinciding with the spectrum of the initial Vq except for the point 
E = oíi = Ek which is absent. At the energy E = a 2 = E\ except for an eigenfunction 
of h\ there exists an associated function (see e.g. 0] and also jlj) which we will also call 
"background eigenfunction". It satisfies the inhomogeneous equation 

Oi - Ei)xe 1 = (Pe, XE t (a) = Xmib) = (10) 

and also the homogeneous one with the squared Hamiltonian 

Oi - Ei) 2 XE t = Xe 1 {o) = XE l (b) = xe 1 + <fE t ■ (11) 

We would like to stress that the set {<^ n }, n — 0, 1, 2, . . .; n ^ k (Ek is "deleted") is not complete 
in L 2 (a, b). To have a complete set one has to add to this set the function XE t 0E] 

As it was already pointed out all spectral points E n , n ^ k, l of ho are spectral points 
of h\ also. So, to prové our claim it remains to analyze only the points E = a\ = Ek and 
E = a 2 = El 

One of the solutions (p^ = (p ai of the Schrodinger equation with E = a± = Ek is given by 

(JüJ) from which it follows that (p^ (a) ^ and (p^(b) ^ 0. A solution vanishing at one of the 
bounds, for instance at x = a 




does not vanish at the other bound. This means that E = a\ = Ek is not a spectral point of 
h x . 

To get a solution of the Schrodinger equation at E = a 2 = E t one can use formula (JÜJ) with 
ipE — ^Ei which gives us the function íçe 1 satisfying the zero boundary conditions meaning as 
it was already mentionned that E — Ei is the spectral point for h\. Moreover, since Ui(a) = 0, 
from (jüj) one gets Pe{.o) = (E — a 2 )ipE(o)- Now if ipE(a) Í>Ei{a) when E — ► Ei, remembering 
that iPe{o) is an analytic function of E having a simple zero at E = E\ (see e.g. [5J we 
conclude that the function PeÍp) is also an analytic function of E but it has a double zero at 
E = a 2 = E[. In such a case together with the function (fE t there exists an associated function 
XEi = (9(Pe/ 9e)e=Ei (see e.g. jSJE])- It is evident that XeX ) = XEiib) = and the equation 
(fTÜ|) it satisfies can be obtained by taking the derivative of equation (0) with respect to E. 
Since íçe = LapE an d L is independent of E one has XEi = Lip El , 4 ! e 1 = {P^eIQe)e=e v The 
function ipEi satisfies the equation (ho — Ei)ipE l = ^e % but it does not satisfy the zero boundary 
conditions which agrees with the fact that ho is a diagonalizable Hamiltonian. Operator L 
(JüJ) turns ipE l into a solution of the equation (fTüj) satisfying the zero boundary conditions thus 
transforming it into a "background eigenfunction" of h\. 

In contrast to the usual SUSY scheme the opposite process, the "deletion" of the level E = E\ 
does not actually delete this level. If we take the Hamiltonian h\ as the initial Hamiltonian 
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for the next second order transformation, leading to the Hamiltonian h 2 , and choose one of 
the transformation functions defming the transformation operator L^ 2 ' of the next step to be 
equal to <pe v actual eigenfunction at E = E\ is deleted but the associated function xe l "comes 
out of the background" and becomes a true eigenfunction of h 2 at E — E\. This statement is 
readily seen if one acts by L^ 2 \ which is constructed in a similar way as L = L^ 1 ' given in 
and intertwines now h\ and h 2 , on both sides of Eq. fjlOjl . takes into account the intertwining 
relation L^ 2 >h\ = h 2 L^ and the property L&'ipfy = 0. If a non-diagonalizable Hamiltonian has 
only one associated function it is transformed in this way into a diagonalizable Hamiltonian. 

3. The simplest example illustrating the possibilities described above is the boundary value 
problem with the zero initial potential Vq(x) = 0. We will choose a = — n and b = n. 
The solutions of the boundary value problem (Q), (JJJ) is well-known, for instance, its discrete 
spectrum is E = E n = \n 2 , n — 1, 2, — 

Let us choose U\ = sin(Ac) and u 2 = exp(-iBx), A,B G R. Formula (jïïj) gives us the 
following "PT-symmetric Hamiltonian: 

Vx ._ 2A\A>-B>) 

[ cos(Ax) — iB sin(Ar)] 2 

For A — 1 the function U\ coincides with the first excited state of h and for B ^ n/2 this 
potential is diagonalizable with the spectrum E = E n = \n 2 , n = 1,3,4, 5,... and E a2 = B 2 . 
For B = 2 the function w 2 is a complex linear combination of the fourth excited state and 
another solution of Eq. (JTJ with Vq(x) = at the same energy and the level E a2 merges with 
the existing level E = 4 which "goes to the background". The potential (fTïïj) becomes non- 
diagonalizable with the discrete spectrum E = E n = ^n 2 , n = 1, 3, 4, 5, . . . studied in detail in 

® 

Now we would like to illustrate the possibility to transform the non-diagonalizable potential 

(|13j) at A = 1 and B = 2 into a digonalizable one. We choose V\ as the initial potential and 

take u\ = f4 and u 2 = fieft where fi e ft is such that ípi e f t (—Tr) = 0. This yields us the following 
potential: 

(k 2 — l)k 2 — 1 — k 2 cos(2x) + cos(2/íx + 2 un)] , , t , 

y 2 = t ^—n vT 1 k + 1 14 

[K COS{KX + K7T) SlllX — Sm{KX + KTT) cosxj^ 

where we denoted a 2 = k 2 - It is regular Vx G (— vr, 7r) provided 0.5 < k < 1.5, k ^ 1 and has 
the spectrum E — E n — n = 1, 3, 4 . . . and E = a 2 = k 2 . For k — 1 (puft = Lip le f t = 
and to realize the transformation with a 2 = 1 one has to use the solution obtained with the 
help of formula (jüj). This corresponds to the backward transformation from V\ to Vq — 
and hence one gets V 2 = 0. For a real k the potential ((Hj) is real and corresponds to the 
Hermitian (hence diagonalizable) Hamiltonian h 2 = —d 2 + V 2 . So, we have transformed 
the non-diagonalizable Hamiltonian h\ to the diagonalizable h 2 . One can also transform 
h\ into a non-Hermitian diagonalizable h 2 by choosing a complex linear combination of two 
linearly independent solutions of equation (J2J corresponding to the same value of E = a 2 as 
transformation function u 2 . 

Our last example is related with the possibility to enlarge the root subspace corresponding 
to E = 4 of the potential (|13|) at A = 1 and B = 2 from the dimension two till the dimension 
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three. For this aim we take U\ = íçe 1 and u 2 = [9 — exp(— 2ïar)]/[l — 3exp(2ú;)] which yields 
the potential 

_ 25e ix + 324è 2 " + 1350e 3íx + 2500è 4 " + 2025e 5íx 
2 ~ (3 + 25e" + 81e 2ix + 75e 3 ") 2 ' 

2 

It has the spectrum E — ^V, n = 3, 4, 5, — 

We hope that the possibility to transform non-diagonalizable PT-symmetric Hamütonians 
to diagonalizable forms may find application in complex quantum mechanics which is currently 
under developments. 
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